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Abstract 

A microscopic state-of-the-art calculation of the nuclear matrix element for neutrinoless double 
beta decay of 150 Nd with an account for nuclear deformation is performed. The proton-neutron 
quasiparticle random phase approximation (QRPA) with a realistic residual interaction [the Brueck- 
ner G matrix derived from the charge-depending Bonn (Bonn-CD) nucleon-nucleon potential] is 
used as the underlying nuclear structure model. The present calculated matrix element is sup- 
pressed by about 40% as compared with our previous QRPA result for 150 Nd obtained with neglect 
of deformation. By making use of this newest nuclear matrix element, one may conclude that neu- 
trinoless double beta decay of 150 Nd, to be measured soon by the SNO+ collaboration, provides 
one of the best probes of the Majorana neutrino mass. 
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Neutrinoless double beta decay (OvfiB decay) is a second order nuclear weak decay process 
with the emission of two electrons only [l|-|3j: (A, Z) — > (A, Z+2)+2e~ . This process violates 
the total lepton number conservation and is therefore forbidden in the standard model (SM) 
of electroweak interaction. The existence of Oz//3/3 decay will immediately prove the neutrino 
to be a Majorana particle (i.e., identical to its antiparticle). 

Furthermore, a study of decay is an indispensable mean to probe the absolute 

neutrino masses at the level of tens of meV. The fact that the neutrinos are massive particles 
was firmly established by neutrino oscillation experiments, thus providing the first evidence 
for physics beyond the SM (for reviews see, e.g., Ref. Q). However, the observed oscillations 
cannot in principle pin down the absolute scale of the neutrino masses. This calls for 
alternative ways one of which is Of/33 decay. 

Thus unambiguous observation of 0v33 decay would be of paramount importance for 
our understanding of particle physics beyond the SM. The next generation of O^/3/3-decay 
experiments (CUORE, GERDA, MAJORANA, SNO+, SuperNEMO, and so on, see, e.g., 
Ref. j3[ for a recent review) has a great discovery potential. Provided the corresponding 
decay rates are accurately measured, knowledge of the relevant nuclear matrix elements 
(NME) M 0v will become indispensable to reliably deduce the effective Majorana mass from 
half-lives Ty 2 of the decay. 

One of the best candidate nuclei for searching decay is 150 Nd since it has the second 
highest endpoint, =3.37 MeV, and the largest phase-space factor for the decay (about 
33 times larger than that for 76 Ge, see, e.g. [l|). The SNO+ experiment at the Sudbury 
Neutrino Observatory will use an Nd-loaded scintillator to search for 0v3f3 decay by looking 
for a distortion in the energy spectrum of decays at the endpoint jBJ. SNO+ will be filled 
with 780 tons of liquid scintillator. The planned loading of 0.1% of the natural Nd translates 
into 43.6 kg of the isotope 150 Nd. It is expected to achieve the sensitivity of T®Y 2 ~ 5 • 10 24 
years after one year of running, with the best final value of about three to four times longer 
(without enrichment of the dissolved Nd). 

Now, to translate the anticipated experimental sensitivity to the decay rate into the sen- 
sitivity expected for the effective Majorana neutrino mass mpp, one needs the corresponding 
NME M 0u . With the result M 0u = 4.74 of Ref. @, already the initial phase of SNO+ will 
be able to probe mpp ~ 100 meV, and will finally be able to achieve sensitivity of mpp ~ 50 
meV corresponding to the case of the so-called inverse hierarchy (IH) of the neutrino mass 
spectrum. 

However, 150 Nd is well-known to be strongly deformed, which strongly hinders a reliable 
theoretical evaluation of the corresponding 0^/3/9-decay NME. For instance, it does not seem 
feasible in the near future to reliably treat this nucleus within the large-scale nuclear shell 
model (LSSM), see, e.g., Ref. Q. Also, the "optimistic" NME of Ref. @ was obtained 
within a microscopic approach, the proton-neutron quasiparticle random phase approxima- 
tion (QRPA), with neglect of deformation. 

Recently, more phenomenological approaches like the pseudo-SU(3) model [gj], the Pro- 
jected Hartree-Fock-Bogoliubov (PHFB) approach |{| and the interacting boson model 
(IBM-2) [To| have been employed to calculate M 0u for strongly deformed heavy nuclei (a 



Ref. [ill]). The results of these models generally reveal a substantial suppression of M 0l/ for 
150 Nd as compared with the QRPA result of Ref. J6J where 150 Nd and 150 Sm were treated as 
spherical nuclei. The recent result of the PHFB [9| is in a fair agreement with the pseudo- 
SU(3) one of Ref. @], but they both are about 1.5 times smaller than M 0u of the IBM-2 [Io|. 
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These results for M 0u give a factor of 2-3 worse limits (as compared with the result of 
Ref. 0]) on the Majorana neutrino mass to be achieved at SNO+ , and basically leave no 
hope to probe the IH region by the current configuration of SNO+. 

Such a spread in calculated NME M 0u for 150 Nd makes it very important to have a 
reliable estimate of the effect of nuclear deformation on M 0u . The most microscopic way 



up-to-date to describe this effect in 150 Nd and 150 Sm is provided by the QRPA. In Refs. 12 
14| a QRPA approach for calculating the 2z//3/3-decay NME M 2v in deformed nuclei was 



developed. The 2z^/3/3-decay half-lives have accurately been measured for a dozen nuclei 
and the corresponding nuclear matrix elements M 2 £ p were extracted (l5j . A theoretical 
interpretation of these matrix elements provides a check of the reliability of different models. 



It was demonstrated in Refs. 12Ml4| that deformation introduces a mechanism of suppression 
of the M 2v matrix element which gets stronger when deformations of the initial and final 
nuclei differ from each other. A similar dependence of the suppression of both M 2u and 
M 0v matrix elements on the difference in deformations was found in the PHFB [9] and the 
LSSM Q. 

In this Rapid Communication we report on the most microscopic state-of-the-art calcula- 
tion of M 0u for 150 Nd with an account for nuclear deformation. The QRPA with a realistic 
residual interaction (the Brueckner G-matrix derived from the Bonn-CD nucleon-nucleon 
potential) is used. The present calculation shows a suppression of M 0u by about 40% 
as compared with our previous QRPA result [6| for 150 Nd that was obtained with neglect of 
deformation. Making use of this newest NME, one may conclude that 0z^/3/? decay of 150 Nd, 
to be searched for by the SNO+ collaboration soon, provides one of the best sensitivities to 
the Majorana neutrino mass and may approach the IH region of the neutrino mass spectrum. 

The NME M 0l/ for strongly deformed, axially symmetric nuclei can be most conveniently 
calculated within the QRPA in the intrinsic coordinate system associated with the rotating 
nucleus. This employs the adiabatic Bohr-Mottelson approximation that is well justified for 
150 Nd, which indeed reveals strong deformation. As for 150 Sm, the enhanced quadrupole 
moment of this nucleus is an indication for its static deformation. Nevertheless, the exper- 
imental level scheme of 150 Sm does not reveal a clear ground-state rotational band. In this 
work we treat 150 Sm in the same manner as 150 Nd. However, a more elaborated theoretical 
treatment going beyond the simple adiabatic approximation might be needed in the future 
to describe the nuclear dynamics of this nucleus. 

Nuclear excitations in the intrinsic system \K V ) are characterized by the projection K of 
the total angular momentum onto the nuclear symmetry axis (the only projection that is 
conserved in strongly deformed nuclei) and the parity n. In Ref. [jjj the structure of the 
intermediate |0 + ) and |1 + ) states was obtained within the QRPA to calculate 2v(3f3-dec&y 
NME M 2v . Here, the approach of Ref. [jjj is straightforwardly extended to calculate all 
possible IK 7 ") states needed to construct M 0u . 

The matrix element M 0l/ is given within the QRPA in the intrinsic system by a sum of 
the partial amplitudes of transitions via the intermediate states K n 

M 0u = M^iK 77 ) ; M 0u (K n ) = s^O^K*). (1) 

Here we use the notation of Appendix B in Ref. [l6[, a stands for the set of four single- 
particle indices {p,p' ,n,n'}, and O a (K n ) is a two-nucleon transition amplitude via all the 
K n states in the intrinsic frame 

O a (K*) = ]T (0+| C Jc n |^m / )(^m / |^m,)(^m,|4 /Cn ,|0+). (2) 

m,i,mf 
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The two sets of intermediate nuclear states generated from the initial and final ground 
states (labeled by rrii and m/, respectively) do not come out identically within the QRPA. 
A standard way to tackle this problem is to introduce the overlap factor of these states 
(K n m,f\K n mi) in Eq. 02]). Two-body matrix elements s^ e ^ of the neutrino potential in 
Eq. (pQ) in a deformed Woods-Saxon single-particle basis are decomposed over the the spher- 
ical harmonic oscillator ones according to the way described in Ref. (lij : 

(def) = V- \p F JK jpJK (sph) I T\ fo\ 

J pp'nn' Z_^i Z-^l PVp nr )n P'Vp'n'Vn' VpVp'VnVn' ^ '> \ ' 

J VpVp' 
VnVn' 




s [ ;A>{J) = Y,{-\)^ +J+J J \ t t J r \ (p(l),p'(2);J\\O i (l,2)\\n(l),n\2);J) (A) 



where J = \J1J + 1, and 0^(1, 2) is the neutrino potential as a function of the coordinates 
of two particles, with i labeling its Fermi (F), Gamow- Teller (GT), and Tensor (T) parts. 
The particle- hole transformation coefficient F^„ n = El B™ (— l) J "~ Qn C/o , from the 
deformed basis into the spherical harmonic oscillator one is constructed from the single- 
particle decomposition coefficients and B™ n (see Ref. [l4| for details), Cj£q Jn _Q n is the 
Clebsch-Gordan coefficient. 

The particle-hole transition amplitudes in Eq. ([2]) can be represented in terms of the 
QRPA forward and backward amplitudes along with the coefficients of the Bogoli- 
ubov transformation u T and v T [14j : 



{Oj\clc n \K w m f ) = v p u n X™t K * + u p v n Y™t K *, 
(K rrii\c p c n \0f) = u p v n X pr ^ K7T + VpU n Y pn l K7r . 

The overlap factor in Eq. (J2J) can be written as: 

(K*m f \K* mi ) = Y\XT/ K ,X^ ~ Y^Y^jK^iBCSflBCSi) (5) 
lilf 

Representations for Hi^ and the overlap factor (BCSf\BCSi) between the initial and final 
BCS vacua are given in Ref. 12| . 



For a numerical computation of the O^/3/9-decay NME M 0v for the process 
150 Nd— )- 150 Sm+2e~, we have straightforwardly extended the method of Ref. 14J. 



The single-particle Schrodinger equation with the Hamiltonian of a deformed Woods- 
Saxon mean field is solved on the basis of a axially-deformed harmonic oscillator. The 
parametrization of the mean field is adopted from the spherical calculations of Refs. 0, 16 



17J. We use here the single-particle deformed basis corresponding in the spherical limit 
to full (4-6)7kJ shells. Decomposition of the deformed single-particle wave functions is 
performed over the spherical harmonic oscillator states within the seven major shells. Only 
quadrupole deformation is taken into account in the calculation. The geometrical quadrupole 
deformation parameter of the deformed Woods-Saxon mean field is obtained by fitting the 
experimental deformation parameter (3 = yf where r c is the charge rms radius and Q p 
is the empirical intrinsic quadrupole moment. The latter can be derived from the laboratory 
quadrupole moments measured by the Coulomb excitation reorientation technique, or from 
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the corresponding B(E2) values [18[. We take in this work the experimental values (3 = 0.29 
and (3 = 0.19 for 150 Nd and 150 Sm, respectively, which are extracted from the B(E2) values 
as being more accurate. The fitted values of the parameter /3 2 of the deformed Woods- 
Saxon mean field, which allow us to reproduce the experimental /3, are listed in Table [H 
The spherical limit, i.e. fi 2 = 0, is considered as well, to compare with the earlier results of 
Ref. [6|. The procedure of fitting (3 2 adopted here is more consistent than the approximate 
ansatz (3 2 = ft used in Ref. [141 ] . 

As in Refs. 0, 14], 16, [l7j], the nuclear Brueckner G matrix, obtained by a solution of the 



Bethe-Goldstone equation with the Bonn-CD one boson exchange nucleon-nucleon potential, 
is used as a residual two-body interaction. First, the BCS equations are solved to obtain 
the Bogoliubov coefficients, gap parameter and the chemical potentials. To solve the QRPA 
equations, one has to fix the particle-hole g p h and particle-particle q pp renormalization factors 
of the residual interaction (see Ref. [lij for details). As in Ref. [HI, we determine a value of 
g p h by fitting the experimental position of the Gamow- Teller giant resonance (GTR) in the 
intermediate nucleus. Since there is no experimental information on the GTR energy for 
150 Nd, we use for this nucleus the same g P h = 0.90 as fitted for 76 Ge (this value is slightly 
different from the fitted g p h = 1.15 of Ref. [l4J because of a different parametrization of the 
mean field used here). The parameter g pp can be determined by fitting the experimental value 
of the 2^/3/3-decay NME M"q T = 0.07 MeV -1 15j. The unquenched axial- vector coupling 



constant gA = 1-25 is used here. The fitted values of g pp are listed Table [D Note, that the 
more realistic procedure of fitting (3? adopted here also gives us more realistic g pp ~ 1 values 



as compared with those of Ref. [14 



TABLE I: The values of the deformation parameter f3 2 of Woods-Saxon mean field for initial 
(final) nuclei fitted in the calculation to reproduce the experimental quadrupole moment. Also the 
fitted values of the particle-particle strength parameter g pp are listed (the particle-hole strength 
parameter is g p h = 0.90). The BCS overlap factor {BC S f\BC Si) © between the initial and final 
BCS vacua is given in the last column. 



initial (final) 


P2 


9pp 


{BCSi\BSC f ) 


nucleus 








150 Nd (150 gm ) 


0.240 (0.153) 


1.05 


0.52 




0.0 (0.0) 


1.01 


0.85 



Having solved the QRPA equations, the two-nucleon transition amplitudes ([2]) are calcu- 
lated and, by combining them with the two-body matrix elements of the neutrino potential, 
the total NME M 0l/ ([1]) is formed. The present computation is rather time consum- 

ing since numerous programming loops are needed to calculate the decompositions of the 
deformed two-body matrix elements over the spherical ones. Therefore, to speed up the 
calculations the mean energy of 7 MeV of the intermediate states is used in the neutrino 
propagator. Following Refs. 0, [lfi], 17], in this first application of the approach the effects 
of the finite nucleon size and higher-order weak currents are included. Recently, it was 
shown [ItJ that a modern self-consistent treatment of the two-nucleon short-range correla- 
tions (s.r.c.) leads to a change in the NME M 0u only by a few percents, much less than 
the traditional Jastrow-type representation of the s.r.c. does. Therefore, we postpone the 
analysis of the anticipated small effects of the s.r.c. to a forthcoming detailed publication. 
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TABLE II: The matrix elements M 0u for the Ov/3/3 decay 150 Nd-> 150 Sm calculated in different 
models. The final result of this work obtained with account of deformation is given in column 4. A 
result with neglect of deformation is also listed (column 3) for comparison with the earlier result of 
Ref. [6] (column 2). The corresponding half-lives T®Y 2 (in years) for an assumed effective Majorana 
neutrino mass {man) = 50 meV are also shown. 





QRPA [6] 


a this work (/?2 


= 0) b this work 


pseudo-SU(3) [8] PHFB [9] IBM-2 [10] 




5.17 


5.78 


3.16 


1.57 1.61 2.32 


Iffc, 10 25 y 


1.72 


1.38 


4.60 


18.7 17.7 8.54 


((7713/3) = 50 meV) 











"using spherical harmonic oscillator wave functions, no deformation allowed. The radius parameter ro = 1.2 
fm is used here, instead of ro = 1.1 fin of Ref. jfjj 
^using Woods-Saxon wave functions, no deformation allowed. 



In Table HTl the presently calculated NME M 0l/ for 150 Nd is listed (column 4) and compared 
with the calculation results by other approaches. One can see that the NME M 0u of this 
work calculated with neglect of deformation (column 3) agrees well with the previous one 
of the spherical QRPA [6j. A small difference can have its origin in the somewhat different 
approximations involved (use of the Woods-Saxon single particle wave functions and the 
BCS overlap factor, neglect of the s.r.c. in the present work). By including deformation 
(column 4), one gets about 1.8 times smaller NME M 0u . The main origin of the suppression 
can be attributed to a smaller BCS overlap factor in the latter case, that is due to a marked 
difference in deformations between 150 Nd and 150 Sm nuclei (see Table H]). 

Our present NME M 0u for 150 Nd, obtained within the state-of-the-art QRPA approach 
that accounts for nuclear deformation, though smaller than the earlier one of Ref. [6] , still is 
significantly larger than the NME of other approaches (columns 5,6, and 7 of Table ILT|) . The 
0^/?/3-decay half-life T®J 2 corresponding to the Majorana neutrino mass {mpp) = 50 meV is 
more than two times shorter as compared with the most optimistic prediction of the IBM-2 
among the other approaches 1 . It allows us to hope that the SNO+ experiment will still be 
able to approach the inverse hierarchy of the neutrino mass spectrum. 

To conclude, in this Rapid Communication the most microscopic state-of-the-art calcu- 
lation of the nuclear matrix element for neutrinoless double beta decay of 150 Nd with an 
account for nuclear deformation is performed. The proton-neutron QRPA with a realistic 
residual interaction (the Brueckner G matrix derived from the Bonn-CD nucleon-nucleon po- 
tential) is used as the underlying nuclear structure model. The 0v(3f3 decay matrix elements 
M 0u calculated in this work shows suppression by about 40% with respect to our previous 
QRPA result for 150 Nd obtained with neglect of deformation. Making use of this newest 
nuclear matrix element, one may conclude that neutrinoless double beta decay of 150 Nd, to 
be measured soon by the SNO+ collaboration, provides one of the best sensitivities for the 
Majorana neutrino mass. 

The authors acknowledge the support of the Deutsche Forschungsgemeinschaft under 



1 Note that by neglecting the Jastrow-type s.r.c. the IBM-2 NME will get about 20% larger and be in 
rather good agreement with our present result. 
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both SFB TR27 "Neutrinos and Beyond" and Graduiertenkolleg GRK683. 
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